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Abstract 

We consider Backward Stochastic Differential Equations (BSDEs) with generators that 
grow quadratically in the control variable. In a more abstract setting, we first allow both the 
terminal condition and the generator to depend on a vector parameter x. We give sufficient 
conditions for the solution pair of the BSDE to be differentiable in x. These results can be 
applied to systems of forward-backward SDE. If the terminal condition of the BSDE is given 
by a sufficiently smooth function of the terminal value of a forward SDE, then its solution pair 
is differentiable with respect to the initial vector of the forward equation. Finally we prove 
sufficient conditions for solutions of quadratic BSDEs to be differentiable in the variational 
sense (Malliavin differentiable). 
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Introduction 

Problems of stochastic control treated by the crucial tool of backward stochastic differential 
equations (BSDEs) have been encountered in many areas of application of mathematics in recent 
years. A particularly important area is focused around optimal hedging problems for contingent 
claims in models of financial markets. Recently, a special class of hedging problems in incomplete 
financial markets has been considered in the area where finance and insurance concepts meet. 
At this interface problems of securitization arise, i.e. insurance risk is transferred to capital 
markets. One particularly interesting risk source is given by climate or environmental hazards 
affecting insurance companies or big branches of the economy that depend on weather such as 
agriculture and fishing, transportation and tourism. The public awareness of climate hazards 
such as floods or hurricanes is continually increasing with the intensity of the discussion about 
irreversible changes due to human impact. 

BSDEs typically appear in the following setting. On a financial market some small investors 
are subject to an external risk source described for instance by weather or climate influences. 
There may also be big investors such as re-insurance companies that depend in a possibly 
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different way on the same risk source. In this situation market incompleteness stems from the 
external risk not hedgeable by the market assets. One may complete the market either by making 
the external risk tradable through the introduction of an insurance asset traded among small 
agents, or by introducing a risk bond issued by a big agent. In this setting, treating the utility 
maximization problem for the agents under an equilibrium condition describing basically market 
clearing for the additional assets, leads to the determination of the market price of external risk 
through a BSDE which in case of exponential utility turns out to be quadratic in the control 
variable (see [HM06J, [CHIM05] and |CIM04j ). Alternatively, instead of maximizing utility with 
respect to exponential utility functions we might minimize risk measured by the entropic risk 
measure. In this setting we again encounter a BSDE with quadratic nonlinearity, of the type 



where W is a finite-dimensional Wiener process of the same dimension as the control process 
Z, with a generator / that depends at most quadratically on Z, and a bounded terminal condition 
£. In the meantime, the big number of papers published on general BSDEs is rivalled by the 
number of papers on BSDEs of this type of nonlinearity. For a more complete list of references 
see [CSTV05] or |KobOO| . In particular, there are papers in which the boundedness condition on 
£ is relaxed to an exponential integrability assumption, or where the stochastic integral process 
of Z is supposed to be a BMO martingale. 

In a particularly interesting case the terminal variable £ is given by a function g{X?f) at 
terminal time T of the solution process X of a forward SDE 



with initial vector Similarly, the driver / may depend on the diffusion dynamics of X x . 

Via the famous link given by the generalized Feynman-Kac formula, systems as the above of 
forward-backward stochastic differential equations are seen to yield a stochastic access to solve 
nonlinear PDE in the viscosity sense, see [KobOOj. 

In this context, questions related to the regularity of the solutions (X x , Y x , Z x ) of the 
stochastic forward-backward system in the classical sense with respect to the initial vector x or 
in the sense of the stochastic calculus of variations (Malliavin calculus) are frequently encoun- 
tered. Equally, from a more analytic point of view also questions of smoothness of the viscosity 
solutions of the PDE associated via the Feynman-Kac link are seen to be very relevant. 

For instance, Horst and Miiller (see |HM06j ) ask for existence, uniqueness and regularity of a 
global classical solution of our PDE from the analytic point of view. Not attempting a systematic 
approach of the problem, they use the natural access of the problem by asking for smoothness 
of the solutions of the stochastic system in terms of the stochastic calculus of variations. But 
subsequently they work under the restrictive condition that the solutions of the BSDE have 
bounded variational derivatives, which is guaranteed only under very restrictive assumptions on 
the coefficients. 

The question of smoothness of the stochastic solutions in the parameter x arises for instance 
in an approach of cross hedging of environmental risks in {AIP05] . Here the setting is roughly 
the one of an incomplete market generated by a number of big and small agents subject to an 
external (e.g. climate related) risk source, and able to invest in a given capital market. The 
risk exposure of different types of agents may be negatively correlated, so that typically one 
type profits from the risky event, while at the same time the other type suffers. Therefore the 
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concept of hedging one type's risk by transferring it to the agents of the other type in a cross 
hedging context makes sense. Mathematically, in the same way as described above, it leads to 
a BSDE of the quadratic type, the solution (Y x , Z x ) of which depends on the initial vector x 
of a forward equation with solution X x . Under certain assumptions, the cross-hedging strategy 
can be explicitly given in a formula depending crucially on x, and in which the sensitivity with 
respect to x describes interesting quality properties of the strategy. 

In this paper, we tackle regularity properties of the solutions (Y X ,Z X ) of BSDEs of the 
quadratic type such as the two previously sketched in a systematic and thorough way. Firstly, 
the particular dependence on the starting vector x of the forward component of a forward- 
backward system will be generalized to the setting of a terminal condition depending in 
a smooth way to be specified on some vector x in a certain Euclidean state space. We both 
consider the smoothness with respect to x in the classical sense, as well as the smoothness in 
the sense of Malliavin's calculus. 

The common pattern of reasoning in order to tackle smoothness properties of any kind 
starts with a priori estimates for difference and differential quotients, or for infinite dimensional 
gradients in the sense of variational calculus. In the estimates, these quantities are related to 
corresponding difference and differential quotients or Malliavin gradients of the terminal variable 
and the driver. To obtain the a priori estimates, we make use to changes of probability of the 
Girsanov type, by which essentially nonlinear parts of the driver are eliminated. Since terminal 
conditions in our treatment are usually bounded, the exponential densities in these measure 
changes are related to BMO martingales. Known results about the inverse Holder inequality 
allow to show that as a consequence the exponential densities are r-integrable for some r > 1 
related to the BMO norm. This way we are able to reduce integrability properties for the 
quantities to be estimated to a natural level. In a second step, the a priori inequalities are 
used to derive the desired smoothness properties from corresponding properties of driver and 
terminal condition. To the best of our knowledge, only Malliavin differentiability results of this 
type have been obtained so far, with strong conditions on the coefficients restricting generality 
considerably (see [HM06J). 

The paper is organized as follows. In section [1] we fix the notation and recall some process 
properties needed in the proofs of the main body of the paper. Section [2] contains the main 
results on classical differentiability. In sections H] and [5] we give a priori bounds for classes of 
non-linear BSDEs. Section [6] contains the proofs of the theorems stated in Section El Section 
[7| is devoted to the application of the proven results to the forward-backward SDE setting. In 
Section [S]we state and prove the Malliavin differentiability results. 

1 Preliminaries 

Throughout this paper let ($7, J 7 , P) be a complete probability space and W = (Wt)t>o a 
d— dimensional Brownian motion. Let {J~t}t>o denote the natural filtration generated by W, 
augmented by the P— null sets of T . 

Let T > 0, £ be an .^-measurable random variable and / : f2 x [0, T] x R x R d -> R. We 
will consider Backward Stochastic Differential Equations (BSDEs) of the form 



As usual we will call £ the terminal condition and the function / the generator of the BSDE 
([!]). A solution consists of a pair (Y, Z) of adapted processes such that ([!]) is satisfied. To be 




(1) 
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correct we should write // (Z t ,dW t ) or Y%=lSt Z\dW\ instead of J t Z t dW t , since W and Z 
are d— dimensional vectors; but for simplicity we use this notation as it is without ambiguity. 
It is important to know which process spaces the solution of a BSDE belongs to. We therefore 
introduce the following notation for the spaces we will frequently use. Let p G [1, 00]. Then, for 
m G N* 



L p (lR m ) is the space of all progressively measurable processes (Xt)t^[o,T] with values in M m 
such that H-XtllJ,, = E[( J Q T |X s | 2 ds) P/2 ] < 00. 

• 7£ p (M m ) is the space of all measurable processes (Xt)te[o,T) with values in R m such that 

= E[^sup tg [ 0] T] \Xt\) ] < 00 • Note that TZ oc (M. m ) is the space of bounded measur- 
able processes. 

• H p (R m ) is the class of all local martingales X such that ||X||^= E P [(X}|] < 00. 

• LP(R m ;P) is the space of .T^-measurable random variables X : 1 — ► M. m such that ||JC||^ P — 
E P [|X| P ] < 00. We will omit reference to the space or the measure when there is no 
ambiguity. 

Furthermore, we use the notation d t = J^, V = (^-, • • • , ^-) for (t, x) G [0, T] x R d . 
Suppose that the generator satisfies, for a > and b, c > 

\f(t,x,y,z)\ < a (l + b\y\) + ^\z\ 2 . (2) 

Kobylanski has shown in [KobOO] that if £ is bounded and the generator / satisfies ([2|), then 
there exists a solution (Y, Z) G 1Z°° x L 2 . Moreover, it follows from the results in [Mor07j, that 
in this case the process Z is such that the stochastic integral process relative to the Brownian 
motion J Q ZdW is a so-called Bounded Mean Oscillation (BMO) martingale. 

Since the BMO property is crucial for the proofs we present in this paper we recall its 
definition and some of its basic properties. For an overview on BMO martingales see Kaxi) l| . 

Definition 1.1 (BMO). Let M be a uniformly integrable {3~t )-martingale satisfying Mq = 0. 
For 1 < p < 00 set 



\M\\bmo p = sup (E 

T stopping time 



IMoo-MrPLTv 



1/p 



The normed linear space {M : \\M\\bmo p < °°} with norm ||M||ba./0 p is denoted by BMO p . If 
we want to stress the measure P we are referring to we will write BMO(P). 

It can be shown that for any p, q G [1, 00] we have BMO p = BMO, (see |Kaz94j ) . Therefore 
we will often omit the index and simply write BMO for the set of BMO martingales. 

In the following Lemma we state the properties of BMO martingales we will frequently use. 

Lemma 1.2 (Properties of BMO martingales). 

1) Given a BMO martingale M with quadratic variation (M), its stochastic exponential 

£(M) T = exp{M T - ^(M) T } 
has integral 1, and thus the measure defined by dQ = £{M)TdP is a probability measure. 
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2) Let M be a BMO martingale relative to the measure P. Then the process M = M — (M) 
is a BMO martingale relative to the measure Q (see Theorem 3.3 in l r Kaz94l). 

3) For any BMO Martingale, it is always possible to find a p > 1 such that £{M) E L p , 
i.e. if \\M\\bmo2 < - ^(p)> then £{M) E LP (see for example Theorem 3.1 \Kaz9J$ ). Where 

^f(x) = |l + ^ log 2(x-i) }^ ~~ 1 f or a ^ 1 < x < 00 an d verifies lim a ._ +1 + ^f(x) = oo and 



2 Differentiability of quadratic BSDEs in the classical sense 

Suppose that the terminal condition and the generator of a quadratic BSDE depend on the 
Euclidean parameter set W 1 for some n E N* . We will show that the smoothness of the terminal 
condition and the generator is transferred to the solution of the BSDE 

Y t x = - [ T Z*dW s + f T f(s, x, Y s x , Z*)ds, x E K n , (3) 
Jt Jt 

where terminal condition and generator are subject to the following conditions 

(CI) / : Ox [0, T]xK n xMxR d — » R is an adapted measurable function such that f(u>, t, x, y, z) = 
l(u,t,x,y, z) +a\z\ 2 , where l(uj,t,x,y, z) is globally Lipschitz in (y, z) and continuously 
differentiable in (x, y, z); for all r > 1 and (t, y, z) the mapping M. d — > L r , x \—* l(u>, t, x, y, z) 
is differentiable and for all x E M. n 

T \n 

and 



limE P f( / \l(s,x',Y s *,Z x ) - l(s,x,Y s \Z x )\ds 

x'^x l\Jo 



0, 



(C2) the random variables £(x) are Tt~ adapted and for every compact set K C W 1 there exists 
a constant c E R such that svvPxtK ||£(zc)||oo — c '> f° r all p > 1 the mapping W 1 — > L p , 
x i— > is differentiable with derivative V£. 

If (QT]) and (Cj2j) are satisfied, then there exists a unique solution (Y x , Z x ) of Equation ([3]). This 
follows from Theorems 2.3 and 2.6 in [KobOOj. We will establish two differentiability results for 
the pair (Y X ,Z X ) in the variable x. We first consider differentiability of the vector valued map 

x^{Y x ,Z x ) 

with respect to the Banach space topology defined on 7Z P (M}) x L p (M d ). This will be stated in 
Theorem 12.11 A slightly more stringent result will be obtained in the subsequent Theorem 12.21 
Here, we consider pathwise differentiability of the maps 

x^(Y t x (u),Z x (io)) 

in the usual sense, for almost all pairs (u>,t). In both cases, the derivatives will be identified 
with (VY X , VZ X ) solving the BSDE 

VY t x = V£(x)- jfvZ x dW s 

+ jf [d x l(s, x, Y x , Z x ) + d y l{s, x, Y x , Z X )VY X + dj(s, x, Y x , Z X )V Z x + 2aZ x VZ x ] ds. 

(4) 



5 



We emphasize at this place that it is not immediate that this BSDE possesses a solution. In 
fact, without considering it as a component of a system of BSDEs also containing the original 
quadratic one, it can only be seen as a linear BSDE with global, but random (and not bounded) 
Lipschitz constants. 

Theorem 2.1. Assume (CUD and ((&). Then for allp > 1, the function R n -> T^QR 1 ) xL p (R d ) ; 
x i— > (Y x , Z x ), is differentiable, and the derivative is a solution of the BSDE Q). 

Under slightly stronger conditions one can show the existence of a modification of Y x which 
is P-a.s. differentiable as a mapping from R n to R. Let ej = (0, . . . , 1, . . . , 0) be the unit vector 
in R n where the ith component is 1 and all the other components 0. For x £ W 1 and h ^ let 
C(x, h, a) = \[£{x + hei) — For the existence of differentiable modifications we will assume 

that 

(C3) for all p > 1 there exists a constant C > such that for all i G {1, . . . , n}, x, x' € R n and 

M'eR \ {0} 



E 



|£(s + ha) - i{x' + ti ei )\ 2p + |C(x, h, <a) - C(x, ti, ei )\ 2p < C(\x - x'\ 2 + \h - h'\ 2 f. 



Theorem 2.2. Suppose, in addition to the assumptions of Theorem \2.1l that (C3) is satisfied 
and that l(t,x,y,z) and its derivatives are globally Lipschitz continuous in (x,y,z). Then there 
exists a function Q x [0, T] x M n — > (u,t, x) \—> (Y t x , Zf )(lo), such that for almost all lu, 

Y x is continuous in t and continuously differentiable in x, and for all x, (Y t x ,Z x ) is a solution 
ofW- 



3 Moment estimates for linear BSDEs with stochastic Lipschitz 
generators 

By formally deriving a quadratic BSDE with generator satisfying (CHJ) and (C[2]) we obtain a 
linear BSDE with a stochastic Lipschitz continuous generator. The Lipschitz constant depends 
on the second component of the solution of the original BSDE. In order to show differentiabil- 
ity, we start deriving a priori estimates for this type of linear BSDE with stochastic Lipschitz 
continuous generator. For this purpose, we first need to show that the moments of the solution 
can be effectively controlled. Therefore this section is devoted to moment estimates of solutions 
of BSDEs of the form 

U t = C- ^ V s dW s + ^ [l(s, U s , V s ) + H S V S + A s ] ds. (5) 

We will make the following assumptions concerning the drivers: 
(Al) For all p > 1, ( is .^-adapted and we have ( G ^(M 1 ), 

(A2) if is a predictable M. d — valued process, integrable with respect to W, such that f HdW is 
a BMO-martingale, 

(A3) I : n x [0,T] x K x R d -> E is such that for all (u,v), the process l(u,t,u,v) is (T t )- 
predictable and there exists a constant M > such that for all (u),t,u,v), 

\l(u,t,u,v)\ < M(\u\ + \v\), 
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(A4) A is a measurable adapted process such that for all p > 1 we have E[( \A s \ds) p ] < oo. 
Moreover, we assume that (U, V) is a solution of ([5|) satisfying 
(A5) [/ Q T U 2 |V^| 2 ds]2 and J Q T |C/ S 74 s |ds are p-integrable for all p > 1. 

Under the assumptions (Al), (A2), (A3), (A4) and (A5) one obtains the following estimates. 

Theorem 3.1 (Moment estimates). Assume that (A1)-(A5) are satisfied. Let p > 1 and r > 1 
such that £{J HdW)x £ L r (P). Then there exists a constant C > 0, depending only on p, T, 
M and the BMO-norm of f HdW), such that with the conjugate exponent q of r we have 



E 1 



sup \U t \ 2p 



■te[o,T] 
Moreover we have 



V,| 2 ds) < CE P 



ICI 



2pq 2 



+ 



\A s \ds) 



(6) 



W[\ \U a \*ds\+W[\ \V s fds) < CE^ \(\ 2q + 



\AJ 2 ds 



(7) 



The proof is divided into several steps. First let j3 > and observe that by applying Ito's 
formula to e^U 2 we obtain 



e^U 2 



Uf-2 



e? T U 2 
(■r 



eP*U.VAW. 



+ J e Ps [- PU 2 + 2U s {l{s,U s ,V s ) + H S V S + A s ) - \V s \ 2 ]ds. 

By (A2), the auxiliary measure defined by Q = £{H ■ W)t • P is in fact a probability measure. 
Then Wt = Wt — Jq H s ds is a Q-Brownian motion, and 



'T 
i-T 



+ J e^ s [(-(3 + 2M)U 2 + 2M\U S \\V S \ - \V S \ 2 + \U s A s \]ds 



By choosing j3 = M 2 + 2M, we obtain 



e pt U 2 + I e 08 (M\U s \ - \V s \) 2 ds < e pT U% - 2 f e f3s U s V s dW s + [ e Ps \U s A s \ds. (8) 
Jt Jt Jt 

We therefore first prove moment estimates under the measure Q. 

Lemma 3.2. For all p > 1 there exists a constant C, depending only onp, T and M, such that 















E Q 


sup \U t \ 2p 


+ E Q 




< CE Q 


"lCI 2p + 




l te[o,T] 











T 



\A s \ds) 



\2p 



(9) 



Moreover we have 



T 



\UJ 2 ds 



L J0 



+ E U 



WJ 2 ds 



< CE Q \\C\ 2 + / |A s | 2 ds 



(10) 
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Proof. Throughout this proof let C\, C2, ■ ■ ., be constants depending only on p, T and M. 
Inequality © implies 



e&U? < e^Ul -2 e? s 



*U s V s dW s + J e^ s \U s A s \ds, 



(11) 



and (A5) together with the existence of the rth moment for £(f HoIW)t yield Ug\V s \ 2 ds € 
L 1 (Q). Hence, since e" l U 2 is (J-i)-adapted, 



|C| 2 + / e?'\U a A s \ds\F t 



(12) 



Integrating both sides and using Young's inequality, we obtain 

E«[ f U 2 ds] < dE«[|C| 2 + C \U s A s \ds] 
Jo Jo 

< C 1 E Q [Q 2 + 2C 1 [ \A s \ 2 ds} + lE Q [[ U 2 ds], 
Jo 2 Jo 

and hence 

E«[ f T U 2 ds] < C 2 E®[\(\ 2 + f T \A s \ 2 ds\. 
Jo Jo 

Inequality ()12|) . (A5) and Doob's LP inequality imply for p > 1 



E«[sup \U t \ 2p ] < C 3 eQ\(\(\ 2 + [ \U s A s \d S Y 
te[o,T] Jo J 



(13) 



< C A E Q 



\C\ 2p +( sup \U t \ I \A s \ds 
te[o,T] Jo 



By Young's inequality, (sup t6[0jT] W)(/ T \A s \ds) p < ^ sup t6[0jT] \U t \ 2 P + 2C 4 (tf \A s \ds) 2 P, 
and hence 



E«[ sup |^| 2 P] < C 5 E Q |C| 2p + ( / |A s |ds " 
te[o,T] L 

In order to complete the proof, note that (JH]) implies 



(14) 



J e^lVsfds 

<e /3T U%-2 [ T e P s U s V s dW s + 2 / T e /3s M|C/ s ||^|ds+ / T e^|C/ s ||A s |ds. (15) 

Ji 



By Young's inequality, 2 jf e^ s M|C/ s ||F s |ds < \ g e^ s \V s \ 2 ds + 8M 2 eP'Ufds, and hence 



r 2 /? a/T 2 r 



-E c 



T 



Jo 



T 



eP s \V s \ 2 ds < EV T C/| + 8M 2 / e (3s U 2 ds+ e^U 2 + e^ s \A s \ 2 ds 



< C 6 E«[C 2 + / T |^ s | 2 ds] 

JO 



which, combined with (|13p leads to the desired Inequality (jlOp . 

Equation (fT5j) . Young's inequality, Doob's L p -inequality and the Burkholder-Davis-Gundy 
inequality imply 



E ( 



< C 7 E C 



T \P 

e^\VJ 2 ds 



o 



|C| 2p + (T^mve^Ui) + (/ ePU.V.dW.) + (/ r ,v '|F,||.L|d.s 



< C 8 E Q |C| 2p + sup e /3tp |C/ t | 2?, + ( / e^ s U 2 s \V s \ 2 dsy + sup |ff t | p 



te[o,r] 



o 



< C7 8 E° 



|C| 2p + sup e^ tp \U t \ 2p + ( sup 
te[o,T] v te[o,T] 



te[o,T] 



^ 8 \A a \dt 



o 



+ sup |tf t | 2p + ( / e ps \A s \ds 



te[o,T] 

By Young's inequality, 



2y 



sup e pt U t 
te[o,T\ 



e^\V s \ 2 ds) 2 <2C7 8 ( sup e tn Uf \ + 
te[o,T] 



Jt TT 2 



P 1 



2Q 



eP'W.Us 



8 Wo 



which implies 

E Q 



T \P 



< CgE C 



< C 10 E C 



\C\ 2p + sup \U t \ 2p +[ I \A s \ds 

te[o,T] 

T 







2p 



ICI 2p + 



UJds 



2p 



Thus, with Inequality (|14p . the proof is complete. 



□ 



Proof of Theorem 1 3. 11 Notice that by the second statement of Lemma ll.2t the process J HdW = 
f HdW - f H*ds belongs to BMO(Q), and hence -/ HdW also. Moreover, £(f HdW)- 1 = 
£{— J HdW). Consequently, by the third statement of Lemma 11.21 there exists an r > 1 
such that £(H ■ W) T € L r (P) and £(H ■ W)? 1 G L r (Q). Throughout let D = max{||£(# • 
W)T\\h r {P)-> -W)^ 1 ||l p (Q)}> Holder's inequality and Lemma [3T2l imply that for the conjugate 
exponent q of r we have 



E p [ sup \U s \ 2p ] = E Q [£(H-W)^ sup \U s \ 2p ] < DE®[ sup \U S 



2p<Jl 



se[o,T] 



s€[0,T] 



se[o,T] 



< Ci D E c 



|C| 2M +( / |^ s |d S 



2p<? 



= CiDE^f^ •VF)t(|C| 2m + (J T \A s \ds 

< c 2 D-^rE p [\c\ 2pq2 + (y iA s |d 



9 



2 P9 1 



where C\,C2 represent constants depending on p,M,T and the BMO norm of J HdW. Sim- 



ilarly, with another constant C 3 , E p [J Q r \V s \ 2p ds] < C3 L' 1 ^ £ E p [|C| 2p92 + \A s \ds 
and hence ([6]). By applying the same arguments to (fTUl) we finally get ([?]). 



2p<J \ 1 



□ 
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4 A priori estimates for linear BSDEs with stochastic Lipschitz 
constants 



In this section we shall derive a priori estimates for the variation of the linear BSDEs that play 
the role of good candidates for the derivatives of our original BSDE. These will be used to prove 
continuous differentiability of the smoothly parametrized solution in subsequent sections. Let 
((, H, h,A) and (£', H', I2, A') be parameters satisfying the properties (Al), (A2), (A3) and (A4) 
of Section [3] and suppose that l\ and 1% are globally Lipschitz continuous and differentiable in 
(u,v). Let (U, V) resp. (U',V) be solutions of the linear BSDE 



U t = C- [ V s dW s + [ [h(s,U s ,V s )+H s V s + A s ]ds 
Jt Jt 



(16) 



resp. 



ui = c- [ T v;dw s + f 

Jt Jt 



both satisfying property (A5). Throughout let 5Ut 
5A t = At — A' t and Sl(t, u, v) = l\{t, u, v) — ^(t, u, v). 



[l 2 (s,U' s X)+H'X + A' s ]ds 
U t - UL SV t = V t 



V{, 5( = c - C, 



Theorem 4.1 (A priori estimates). Suppose we have for all (3 > 1, Jq 5U 2 \5V s \ 2 ds € L@(P) 
and J Q T \6U s 6A s \ds G L^{P). Let p > 1 and r > 1 such that £(f H'dW) T £ L r (P). Then there 
exists a constant C > 0, depending only on p, T, M and the BMO-norm of J H'dW, such that 
with the conjugate exponent q of r we have 



E p 


' sup |fl7 t |*" 


+ E P 




l t£[0,T] J 





<c{e p [|5C| 2jV + ( J 



\SV s \ 2 ds 



\5l{s,U' s X) + 5A s \ds 



+ (E ; 



\2pq 2 



+ 



T 



A^dsf^'])^^ 1 



\Hc - H'\ 2 ds 



2pq 2 



2,-r 



We proceed in the same spirit as in the preceding section. Before proving Theorem 14.11 
we will show a priori estimates with respect to the auxiliary probability measure Q defined by 
Q = £(J H'dW) T ■ P. Note that W t = W t - J * H' s ds is a Q-Brownian motion. 

Lemma 4.2. Let p > 1. There exists a constant C > 0, depending only on p, T and M, such 
that 





' sup \SU t \ 2 P 


< c{e q 


\S(\ 2p + 




-te[o,T] 










+ (e q 


ICI 2p + 



\2p 



(17) 



,2p 



T 



E c 



[(^ T |^| 2 d s ) P ] < c{mQ[\5C\ 2p +(£ \5l(s,U' s X) + SA s \d. 



\H,-H'J 2 ds 



2p- 



2p- 



(18) 



+ (e«[|C| 2p +(^ L4 s |ds) 2p ]) V[(^ \H s -H' s \ 2 ds) 2P f} 



10 



Proof. The difference SU satisfies 



8( - [ T SV s dW s + [ T [(H S V S - H'X) + h(s, U s , V 8 ) - l 2 (s, U' a , Vj,) + 8A s ]ds 
Jt Jt 



5( 



f T 5V s dW s + [ T [h(s, U' s , V s ') - l 2 (s, U' s , V' s ) + H' S 5V S + 5A s ]ds 
Jt Jt 

+ j\{H s - H' S )V S + h(s, U s , V s ) - h(s, U's, V^ds. 
Let f3 > 0. Applying Ito's formula to e^ l 8U 2 ,t > 0, yields the equation 



e^SU? = e^ T 8Ul - 2 / e Ps 8U s 8V s dW s + 2 / eP'SU. H' 5V s ds 



T 



T 



(55U 2 s ~ Ws\ 2 + 2(/i(s, U a , V s ) - h{s, U's, Vl))5U a 

+2 [ e Ps 5U s (H s -H' s )V s ds + 2 [ eP a 8U a (8l a + 8A.)6s, 
Jt Jt 



ds 



(19) 



where 81 s = h{s, U' s , V s r ) — fo(s, U' s , V a ). Using the Lipschitz property of l\ we obtain 



e^SU? < e pT 8U^+ e Ps \{-p + 2M)8U 2 - \8V S \ 2 + 2M\8U S \ \8V S \ 



ds 



+2 f e^8U s [(H s - H' S )V S + Sis + 8A s ]ds - 2 [ e ps 8U s 8V s dW s 
Jt Jt 



li(3= (M 2 + 2M), then 



e 0t 8U?+ f e^ s (M\8U s \-\8Vs\) 2 ds < e pT 8U^ + 2f e^ s 8U s [(H s - H' S )V S + 8l a + 8A s ]ds 
Jt Jt 



2 I e Ps 8U s 8V s dW s . 



(20) 



We will now derive the desired estimates from Equation (|20p . First observe that by taking 
conditional expectations, we get 



e^8U 2 < e? T ¥P 



8U% + 2^ e^ s 8Us[{H s - H' S )V S + 8l s + 8A s ]ds\F t 



Let p > 1. Then for some constants C\, C%, . . ., depending on p, T and M, we obtain 

sup \8U t \ 2p < Ci sup {(¥Pl\8U T \ 2 \rt] +E[ / \8U S [(H S - H' S )V S + 8l s + 8A s ]\ds\ F t ]Y\ 
te[o,T] te[o,T] I v Jo ' ) 

and by Doob's IP inequality we get 



E Q [ sup \8U t \ 2p ] < C 2 \eP \8U T \ 2p ]+E[[ I \8Us[(H s - H' S )V S + 8l s + 8As]\ds) } 

te[o,T] 
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By using Young's and Holder's inequalities we have 
rT 

\5U S [(H S - H' S )V S + 8l s + 5A s ]\ds 

sup \5U t \ p \( [ T 
te[o,T] LWo 

sup \6U t \ 2p 

t£[0,T] 

■T 



E Q 


U 






< C 3 E Q 






< 


1 E 


" 2C 4 



\H q -H'fds 



\V s \ 2 ds + / \Sl 8 +6A 8 \ds 



+4C 4 E Q 



\H 9 -H'fds 



o 



\V s \ 2 ds)) + 



|5/ s + <L4 s |ds 



2p 



2Ci 



' sup |«5C/*| 2 p 


+ c 5 { 


L te[o,T] J 





2p 



+E' 



Q 



(H s - H' s ) 2 ds 



2p 



iKi 2 d») 2 T} 



(21) 



Therefore, we may further estimate 



E«[ sup TO 2p ] < C 6 {E«[|5C| 2p ]+IE Q [( / |<H fl + <L4 s |ds) 
te[o,T] L Jo 



T 



2 Pl 



-E 1 - 



|# s -ii"ll 2 ds 



2p 



Due to Lemma 



/o |T4| 2 ds) < C 7 E® \(\ 2p + (J |^ s |ds) 2 f < oo, which implies 



the 5U S part of Inequality ([T?P . 

In order to prove the second inequality, note that (1201) also implies 



/ e^l^pds < e l3T 5U% + 2 [ e ps 8U s [{H s - H' S )V S + Sl s + tfAJdi 



T 



+2M / e^ s |<5C/ s | |<5V;|ds-2 / e^^C/^^dWs- 
Equation (j22|) . Doob's L p -inequality and the Burkholder-Davis-Gundy inequality imply 



(22) 



T 



E«[( / |<^| 2 d S ] < C 8 (E« |5C| 2p + / |«/,| 2p ds + E«[( / ^ 2 <5|^| 2 d S 2 ] 



+E Q [( / \5U S [(H S - H' S )V S + 5l s + 5A s ]\ds) } 
Jo 



Consequently, Young's inequality allows to deduce 

E Q [( [ T \5V s \ 2 dsY] < C 9 {e q \\5C\ 2p + [ T \5U s \ 2p ds\ +E Q [ sup \5U t \ 2p ] 
Wo ' I L Jo i te[o,T] 

rT 



+^ Q [{J q \SU S [(H S -H' S )V S + Sl s + <L4J|da) p ] J. 



Finally, (JT7D and jH]) imply 

rT 







E w [( / |<5K| 2 ds)"] < Ci E y |5C| 2p + ( / |«. + <5A»|da 



2p- 



+C 10 EQ[|C| 2p + (/ |^|ds) 2p ]iE 



r<3 



is, - fl"'J 2 ds 



2p 
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and hence the proof is complete. 



□ 



Proof of Theorem \4-l\ This can be deduced from Lemma 14.21 with arguments similar to those 
of Theorem 13.11 We just have to invoke Lemma ll. 21 □ 



5 A priori estimates for quadratic BSDEs 

Consider the two quadratic BSDEs 



and 



Y t = C- [ Z s dW s + [ [l 1 (8,Y a ,Z a ) + aZ?]ds 
Jt Jt 

y{ = c'- [ T Z' s dW s + ! [l 2 (sX,K) + ^(K) 2 ]ds, 
Jt Jt 



(23) 
(24) 



where £ and £' are two bounded J^-measurable random variables, and l\ and I2 are globally 
Lipschitz and differentiable in (y,z). Put now SYt = Yt — Y/, 5Zt = Z% — Z' t , = £ — and 
81 = h — h- The a priori estimates we shall prove next will serve for establishing (moment) 
smoothness of the solution of the quadratic BSDE with respect to a parameter on which the 
terminal variable depends smoothly. Note first that by boundedness of £ and £' we have that 
both J ZdW and J Z'dW are BMO martingales, so that we may again invoke the key Lemma 

Ed 

Theorem 5.1. Suppose that for all (3 > 1 we have J Q T \5l(s,Y s , Z s )\ds € L^{P). Let p > 1 
and choose r > 1 such that £(a(Z s + Z' s ) ■ W)t G L r (P). Then there exists a constant C > 0, 
depending only onp, T, M and the BMO-norm of (a f (Z s +Z' s )dW), such that with the conjugate 
exponent q of r we have 



E 1 



2pq 2 









' sup \SY t \ 2 P 


+ E P 




l te[o,T] 







<C E p 



2pq 2 



+ 



\5l(s,Y s ,Z s )\ds) 



Moreover we have 

cT 



E 1 



\5Y s \ 2 ds] +E J 







Jo 



\SZJ 2 ds] < C E 



m 2q + 



\5l{s,Y s ,Z s )\ds 



2r 



We give only a sketch of the proof since the arguments are very similar to the ones used in 
the proofs in Sections [3] and [U 
First observe that 

5Y t = 5^- f T 5Z s dW s + i' T [l l {s,Y s ,Z s )-h(sX,Z' s ) + SKsXiZ's)+a{Z s + Z' s )5Z s ]ds. 
Jt Jt 

By applying Ito's formula to e^*|<5Yi| 2 we obtain 
e^\5Y t \ 2 -e^ T \5Y T \ 2 

= 2^ e^F s (Zi( S) y s ,Z s ) - h(sX,Z' s ) + Sl{sX,Z' s ))ds -2 J eP s (P\SY a \ 2 + \5Z s \ 2 )ds 

+2 [ e^ s 5Y s a(Z s + Z' s )5Z s ds - 2 [ e ps 5Y s 5Z s dW s . (25) 
Jt Jt 
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We start with a priori estimates under the auxiliary probability measure Q defined by Q 
S(a J (Zg + Z' s )dW) ■ P. Note that W t = W t - J * a(Z s + Z' s )ds is a Q-Brownian motion. 
Let f3 > 0. Equality (|25p and the Lipschitz property of /i yield 

e pt \8Y t \ 2 < e f3T \5t\ 2 -2j e^ s SY S SZ s dW s + 2 J e^ s SY s 5l(s,Y s \Z' s )ds 

+ J e l3s ((-(3 + 2M)\5Y s \ 2 - \5Z s \ 2 + 2M\5Y s \\6Z s \^ds. 

By choosing (3 = M 2 + 2M we obtain the general inequality 

e^\5Y t \ 2 + J T e^ s (\6Z s \ - M\5Y s \) 2 ds 



T 



< e pi + / e ps 5Y s Sl(s, F s , Z' s )ds - 2 / e>" s 5Y s 6Z s dW s . (26) 
Jt Jt 

Note that the process J* * e@ s 5Y S 5Z s dW s is a strict martingale because 5Y S is bounded and 
(<fZ • W) is BMO relative to Q. 

Notice that Equation (|26|) is of similar but simpler form than Equation (|20p . This is because 
the (H s — H' s ) term in (|20l) has been completely absorbed by the Girsanov measure change. As 
a consequence, following the proof of Lemma 14.21 we obtain the following estimates: 

Lemma 5.2. For all p > 1 there exists a constant C > 0, depending only on p, M and T, such 
that 









' sup \5Y t \ 2 P 


+ E Q 




L te[o,T] J 







Moreover we have 

rT 



r r T 

< CE Q m 2p + ( \Sl(s,Y s ,Z t 
1 Jo 



)\ds) 



2p 



E Q [[ \5Y t \ 2 ds] +E Q [ [ \5Z s \ 2 ds] < CE Q 
Jo Jo 



^| 2 + / \d~l(s,Y s ,Z s )\ds 



■ (27) 



Proof of Theorem \5.1l The arguments are similar to those of the proof of Theorem 13.11 Just 
make use of Lemma 11.21 □ 



6 Proof of the differentiability 

We now approach the problem of differentiability of the solutions of a quadratic BSDEs with 
respect to a vector parameter on which the terminal condition depends differentiably. We start 
with the proof of the weaker property of Theorem 12.11 Our line of reasoning will be somewhat 
different from the one used for instance by Kunita |Kun90] in the proof of the diffeomorphism 
property of smooth flows of solutions of stochastic differential equations. He starts with formally 
differentiating the stochastic differential equation, and showing that the resulting equation pos- 
sesses a solution. The latter is then used explicitly in moment estimates for its deviation from 
difference quotients of the original equation. The estimates are then used to prove pathwise 
convergence of the difference quotients to the solution of the differentiated SDE. We emphasize 
that in our proofs, we will have to derive moment estimates for differences of difference quotients 
instead. They will allow us to show the existence of a derivative process in a Cauchy sequence 
type argument using the completeness of underlying vector spaces, which of course will be the 
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solution process of the formally differentiated BSDE. So our procedure contains the statement 
of the existence of a solution of the latter as a by-product of the proof of the Theorem EU It is 
not already available as a good candidate for the derivative process, since, as we stated earlier, 
the formally differentiated BSDE is a globally Lipschitz one with random Lipschitz constants 
for which the classical existence theorems do not immediately apply. Throughout assume that 
f(t,x,y,z) = l(t,x,y,z) + a\z\ 2 and £(x) satisfy (CTT]) and (CO]) respectively. 

For all ieR" let (Yf , Zf) be a solution of the BSDE ©. It is known that the solution is 



unique and that {Y x , Z x ) G TZ°° (R 1 ) x L 2 (R d ) (see [KobOOj ). 

It follows from Lemma 1 in [Mor07] that there exists a constant D > such that for all 
x G R n we have \\(Z X ■ W)t\\bmo 2 — D. Now let r > 1 be such that \&(r) > 2aD (see property 
3) of Lemma I1.2D . and denote as before by q the conjugate exponent of r. 



Proof of Theorem \2.1[ To simplify notation we assume that M > is a constant such that 
x G R n , and the derivatives of I in (y,z) are all bounded by M. We first show that all the 
partial derivatives of Y and Z exist. Let x G R n and e,{ = (0, . . . , 1, . . . 0) be the unit vector in 
R n the ith component of which is 1 and all the others 0. For all h ^ 0, let Up = j.(y t x+eih - Y t x ), 



= l( Z ?+ h «< - Z x ) and C h = + he i) ~ *(*))■ 

Let p > 1. Note that for all h ^ 

U t h = C h -J t T V s h dW s + £ -[/(*, x + he h Y? +h «, Z 



x+hei 



f(s,x,Y x ,Z x )]ds. 



1 



To simplify the last term we use a line integral transformation. For all (a>, t) G f2 x R + let s Xj h = 
s x>h {u,t) : [0,1] -^R"+!+ d be defined by s x>h {9) = (x+0he h Y t x +6(Y x+he * -Y t x ), Zf +9(Z x+fiei - 
Z x )). Though s Xt h depends on i we omit to indicate this dependence for ease of notation. Note 
that y Xth (0) = (e u U t h ,V t h ). Moreover, A^ h = £ £(8 x , h (e))M, G x t ' h = § ( 
It' h = lo §^( s x,h(6))d6 are (^f)-adapted processes satisfying 

l(t,x,Y?,Zf)] = [\vi( Sx , h (9)),s' h (9))de 



-(s x , h (6))de and 



-[l(t,x + hei,Y t 



x+hei yx+hei\ 



7 1 







A x,h + G x, nu n + j. 



~iX,h TT h 



x,h-yh 



Since the derivatives of I are bounded by M, G x,h and I x ' h are bounded by M as well. However, 
we stress that A x,h is not necessarily bounded. We define two random functions m x ' h {u,v) 



x.hi 

m s [u,v) 



(Gs ,h u + Is ,h v) and m s (u,v) 



and m s (u,v) from R + to R such that (it, u) 
[d y l(s, Y x , Z x )u + d z l(s, Yf, Z x )v\. Observe that these functions satisfy (A3) and that they are 
Lipschitz continuous and differentiable in (u,v). In these terms, 



U 



^_ [ T Vs h dW s + r[m x > h {UtV^)+A x > h + a(Z ; 
Jt Jt 



x+hei 



Z x )V s h ]ds, 



and thus we obtain an equation as modelled by ((5|). Notice that for all h, h' ^ the pairs 
(U h ,V h ) and (U h - U h ' ,V h - V h ') satisfy assumptions (A4) and (A5). Therefore Theorem ETT1 



implies with 6 At = A 



x,h 



A 



x,h' 



E[ sup \U t 
L te[o,T] 



h 



h'\2p 



<c{e[\c 
+^[\C h ' 



h 



c 



ti\2pq 2 



+ 



2pq 2 



+ 





x,h 



\Ar |da J 



m x ' h (U ; 

\2pq 2 



E 



m 



x',h' 



(U*,V s h ) + 5A s \ds 



2pq 2 



1 



a 2 \Z 



x+hei 



Z 



x+h'ei |2 



(28) 



*) ]*}. 
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Condition (CED implies that E[\( h - C h '\ 2pq } converges to zero as h, h! — ► 0. Moreover, for some 
open set O containing we have sup h / eC> ^ y(E,\( h \ 2pq ) < oo. Due to Condition (QlD), we may 

also assume that sup ft / gC ,y{ } E[(J T \A x ' h \ds) 2pq2 ] < oo. Moreover, 

r T 

limE( / \l(s,x + hei,Y s x ,Z*)-l(s,x,Y s x ,Z*)\dsf = 
h->o J Q 

for all ft > 1, and therefore, with Theorem 15.11 the third summand on the right hand side of 
([28]) converges to zero as h, b! — > 0. 

In order to prove convergence of the second summand let P®\ be the product measure of P 
and the Lebesgue measure A on [0, T\. It follows from Theorem 1 5 . 1 1 that Z x+hei converges to Z x 
in measure relative to P <8> A. Moreover, for all t S [0, T], Y x+he ' converges to Y x in probability. 
Since the partial derivatives l y and l z are continuous and bounded, dominated convergence 

implies lim.f l ^ E p lf Q \m x ' (Ug,V^)—m s (Ug,V^)\dsj = 0. Condition (Ql]) guarantees 



lim/ l ^ E p ( J \5A s \dsj = 0, and hence, lim/^/^oE sup ie [ T ] \U^ - U^' \ 2p 
Finally, Theorem 14.11 and an estimation similar to (|28p yield 

lim E( / W! 1 - V Q h '\ 2 ds 1 -(). 

h,h'-*0 

Now let (h n ) be a sequence in R\ {0} converging to zero. Then, since lZ 2p (R}) and L 2p (M d ) 
are Banach spaces, the sequence U hn converges to a process J^-Y x , and V h " to a process 
-§^r.Zf with respect to the corresponding norms. By convergence term by term for the difference 
quotient version of the quadratic BSDE and its formal derivative, which follows from our a priori 
estimates, we see that the pair (-£^-Y x , jj^rZf) is a solution of the BSDE 

dY x d w s f T dZf rT 



dx; dx 



f r)7 x f r 

Z(x)-J -g^dW s + J [d x d(s,x,Y x ,Z x ) 

ciY x r)7 x B7 1 

+d y l{s, x, Y x , Z x )-^ + d z l(s, x, Y x , Z x )-± + 2aZ x * 



d Xi * v ' ' s ' 3 ' d Xi s dxi 



ds. 



sup te[0)T ] \U t h - £-Y?\ 2p 



and 



Similarly to the first part one can show that lim^o E 

lim h ^ E^ \V s h - £-Z x \ 2 d s y = 0, and thus R n -» ft^QR 1 ) x h 2p (R d ), x i-> {Y t x ,Zf) is 
partially differentiable. The a priori estimates of Theorem 14.11 imply that the mapping x \— > 
(X7Y X ,X7Z X ) is continuous and hence, (Y t x , Zf) is totally differentiable. Since differentiability 
with respect to 2pth moments implies differentiability with respect to all inferior moments above 
1, we have established the result. □ 

As a byproduct of the previous proof we obtain that for every x £ W 1 there exists a solution 
(VY t x ,VZf) of the BSDE ©. 

We now proceed with the proof of Theorem 12. 21 in which we claim pathwise continuous 
differentiability. To be consistent with the previous proof, we will again compare difference 
quotients varying in h. To this end we need the following estimates. 

Lemma 6.1. Suppose (C3) is satisfied and that I and the derivatives of I are all Lipschitz 
continuous in (x,y,z). Then for all p > 1 there exists a constant C > 0, dependent only on p, 
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T, M and D, such that for all x, x' £ W 1 , h,h' £ R and i £ {1, . . . , n}, 

2p 



E 



sup 

te[o,T] 



x+hei 



Y 



x'+h'ei 



+ 



T \ P 

x+hei ^x'+h'ei |2 



4 ds 



Proof. This follows from Theorem 15.11 where we put l±(s,y,z) = l(s,x + hei,y,z), l2(s,y,z) = 
l(s,x' + h'ei,y,z). □ 

The preceding Lemma immediately implies a first pathwise smoothness result in x for the 
process Y x . In fact, Kolmogorov's continuity criterion applies and yields a modification of Y x 
which is continuous in x. More precisely: 

Corollary 6.2. There exists a process Y x such that for all (t, to) £ [0, T] x 0,, the function 
x i — ► Y x (uj) is continuous, and for all (t,x) we have Y x = Y x almost surely. 



Let e; be a unit vector in R n . For all x £ W 1 and h ^ 0, let U 7 



x.h 



1 / v x+hei V x 



v, 



x,h ^ l^x+het _ and C x,h = + hei) _ ^ x) y j f h = Q; then define ^,0 = _d_y X ^ 

x,0 



dxi 



Z x and C x, ° = jf~£( x )- The proof of Theorem l2.2l will be based on the following result 



on the usual difference of difference quotients. Knowing a "good candidate" for the derivative 
from Theorem 12.11 we allow h = this time, by replacing the difference quotient with this 
candidate. 

Lemma 6.3. Let p > 1 and O C M n+1 be an open set contained in a ball of radius k. Suppose 
that Condition (C3) holds and that I and the derivatives of I in (x,y,z) are Lipschitz continuous 
in (x,y,z) with Lipschitz constant L > 0. Then there exists a constant C, depending on n, L, 
p, T, M , D, such that for all (x, h) and (x' , h') £ O, 



E 



sup \U t x ' h - Uf' h '\ 2p ] < C(\x - x'\ 2 + \h- h'\ 2 f. 



te[o,T] 



(29) 



Proof. Throughout the proof, C\, C2, . . . are constants depending on k, L, p, T, M, D. 

Since O is bounded, (C3) implies that for every r > 1 there exists a constant C\ such that for 
all (x, h) £ O we have E(sup ie[0>T] |CT' h | 2r ) < d- Now let s x>h , m x > h , A x > h , G x > h , I x ' h and U x ' h be 
defined as in the proof of Theorem 12 .2\ and denote A x '° = §^(x, Y x , Z x ), G x '° = f^(x, Y x , Z x ), 
etc. Then the estimate (I29p will be deduced from the inequality 



E 

< C 2 lE 
+ E 

+ E 



sup \U t ' - U t I p 
te[o,r] 



(30) 



x.h 



c 



x',h'\2pq 2 



x '> h '(U x > h , V x ' h ) - m x ' h (U x ' n , V* 



^x,h ujx,h T/x,h\ 



+ \Af h -A x ' n \ds\ 



x'h' 







\2p<? 2 



+ ( I* \A x '> h '\ds) 2pq ^^ 



E 



T 



a 2\ Z x'+h'e l 



z 



x+hei 1 2 



dsj 



2pq 2 



-9- 



which follows from Theorem 14. 1[ We first analyze the order of the convergence of 



Bi(x, x , h, h') = E 



\m x s > h {Uf\ V x,h ) - m x ' h (U x ' h , V x ' h )\ds 



1 

71 



as h, h' — > 0. 
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To this end notice that 

Bi.{x,rf ,h,h') < 



C 3 { U \G*'' h ' -G*' h \\U?> h 



\dt 



2pq 2 \ 



+ j E ( / |if' h '-if' ft [|Vf ,ft |dt 



Then 



f T \G?' h ' -C%' h \\U?' h \dt 
Jo 



< sup I £/ f 

te[o,T] 



\G( > h ' - G*' h \dt 



\d y i{s x , M {o))-d y i{s x de))\de)&t 



< sup r i r 

te[o,T] Jo \Jo 

< sup \u?> h \ f T f 1 L\s x , }h/ (o)-s Xjh (e)\dedt 

telO.T] JO JO 



<C 4 sup \U^ h \(\x'-x\ + \h'-h\+ sup \Yf -Y t x \+ sup \ Y x '+ h '^ - Y x+he ^\ 
te[o,T] ^ te[o,T] te[o,T] 

+ J (\zf -Zf\ + \ zf +h ' e > - Z x+he * \)dt), 

and, by applying Holder's inequality we obtain with Lemma 16, II 



E 



Similarly, ( E 



\G?' h ' - G x ' h \\U?' h \dt 



2pq 2 \ 



< C 5 (\h-h'\ 2 + 



x — X 



>\2\P 



rT I T x'h' T X,h I irrX.h i i. 

Jo \h ' -It Wt 

elude B x (x, x', h, h') < C 7 (\h - h'\ 2 + \x - x'\ 2 ) p 
By using similar arguments we get 



2pq 2 \ ? 



< 6*6 — h'\ 2 + \x — x'| 2 ) p , and so we con- 



E 



T 



\A x ' h -A x '' h '\dt 





1 






< E 











1^(5*^0)) -0^,/.(0))ld0dt 



< C 8 E 




|%^W-^W|d0di] 



2pq 2 



'0 JO 

,'|2 



< c 9 (\h-h'\ 2 + \x-x'\ 2 y. 



Theorem 15.11 and the Lipschitz continuity of I imply 

^ 2pq 2 

e' 



T \ • 

\r?x+hei _ Z x ' +fl ' ei \ 2 ds 



Z: '"' -V/al 

< CiqE [I (e(x + to*) - i{x' + /i'e,) | 4 ^ 

+ (J T \l(s,x + hei, Y x+he \ Z x+he <) -l(s,x' + h' ei , Y x+he \ Z x+he *)\ds} 

< C n (\x-x'\ 2 + \h-h'\ 2 y. 



2q* 
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Finally, (C3) yields (E\( x > h - ,h' < C 12 {\x - x'\ 2 + \h- h'\ 2 )'P, and hence 



E 



sup | U t 

te[o,T] 



X,tl 



x',h'i2p 



< c 13 (\x - x '\ 2 + \h - h'\ 2 y. 



□ 

Proof of Theorem \2.S[ To simplify notation we may assume that (|29p is satisfied for O = R n+1 . 
Assume that Y x is continuous in x (see Corollary 16. 2p . Lemma [6. 3 1 and Kolmogorov's continuity 
criterion imply that U^' has a modification fjf' continuous in (x,h). Define -j^rY® = U^'° 
and note that we obtain thus a continuous version of the solution of the BSDE (J2J). For all 
(x, h) G Q n+1 let N(x, h) be a null set such that for all u (£ N(x, h) we have U?' h (uj) = Vf' h {cj). 
Then, N = U(a;Me<Q> n + 1 N{x,h) is a null set such that for all to ^ N the following implication 
holds: If q k G Q n and r k G Q \ {0} are sequences with lim^oo q k = x G M n and lim^oo = 0, 
then 

lim l(Y t qk+rkSi - Y t qk ) = ^-Y x . 

As a consequence of this and the subsequent Lemma 16.41 Y x (tu) is continuously partially differ- 
entiable relative to x% if u ^ N. Since we can choose such a null set for any i G {1, . . . , n}, total 
differentiability follows and the proof is complete. □ 



Lemma 6.4. Let f : R n — > M 6e a continuous function and g : M n — > R n a continuous vector 
field. Suppose that for all sequences q k € Q n with q k — > x G K™ and r fc G Q \ {0} uirf/i — > 

1 



- /(%)) = ffi(^), 
where 1 < i < n. Then f is differentiable and V/ = g. 



lim —{f{q k +r k ei 
k — >oc r k 



Proof. To simplify notation assume that n = 1. Let i^el with x k —* x G M and /i^ G M \ {0} 
with /ifc — > 0. Since / is continuous we may choose q k G Q and G Q \ {0} such that 
\f(q k ) ~ /(x fc )| < |/(g* + r fc ) - f(x k + /*)] < J^l and 1 1 



2 k 



J_| < J_ Then 



|r- (f{x k + hk) - f{xk)) - g{x)\ 



< 



^[{f(x k + h k ) - f{x k )) - (f(q k + r k ) - f{q k ))] 



+ 



1 

h k 



-)(/(% + r k ) - f(q k ))\ + \-{f(q k + r k ) - f{q k )) - g(x)\ 
r k r k x 



0, 



(A; — * 00), 



and hence / is partially differentiable. Since the partial derivatives gi are continuous, / is also 
totally differentiable. □ 



7 Differentiability of quadratic Forward-Backward SDEs 

In this section we will specify the results obtained in the preceding sections to BSDEs where 
the terminal conditions are determined by a forward SDE driven by the same Brownian motion 
as the BSDE. When considering BSDEs with terminal condition determined by a forward SDE 
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we will need regularity of the forward equation. This will be guaranteed if the coefficients are 
functions belonging to the following space. 

Throughout this section let again n be a positive integer and W a d-dimensional Brownian 
motion. 

Definition 7.1. Let k, m > 1. We denote by B kxm the set of all functions h : [0,T] xl"^ 
R fcxm , (t,x) i— > h(t,x), differentiate in x, for which there exists a constant C > such that 



su P(t,x)e[o,T]x 
x 



dh{t,x) 



< C; and for all t G [0,T] we have sup,^ 



l+\x\ 



< C and 



* dh dx^ ls Lipschitz continuous with Lipschitz constant C. 
With any pair h G B nxl and a G B nxd we associate the second order differential operator 



dxidxj ' 

We will consider Forward-Backward SDEs (FBSDEs) of the form 

J X? = x + £ b(s, X x )ds + /„* a(s, Xf)dW t , x G 1 

9 {x x ) + ; T /( Sj x* y/, z-)d s - jf z x dw s , 



VI 



(31) 



where the coefficients satisfy the following assumptions: 
(Dl) a G B nxd , b G B nxl , 

(D2) / : Q x [0, T] x R™ x R x R d -> R and 5 : R n -> R are measurable functions such 
that f(u,t,x,y,z) = l(u>,t,x,y,z) + a\z\ 2 , where t, x, y, z) is globally Lipschitz and 
continuously differentiable in (x,y,z), 

(D3) g : R n — ► R is a twice differentiable function such that Vg ■ a G B lxd and Cg G B lxl . 

It follows from standard results on SDEs and from Theorem 2.3 in [Kob00| that there exists a 
solution (X x , Y x , Z x ) of Equation (|3ip . As we will show, the results of Section [2] imply that 
(X x , Y x , Z x ) is differentiable in x and that the derivatives (VF, VY X , VZ X ) solve the FBSDE 

VXf = 1 + Jj d x b(X*)VX?d8 + £ ^(X^VXfdW,, 

vyf = d x g(x x )vx x - jf vz x dw s + jf \dj( s , x x ,y x , z x )vx x + d y i( s , x x , y x , z x )vy x 



+d z l(s, X x , Y s x , Z X )VZ X + 2aZ x VZl 



ds. 



(32) 

Our first result parallels Theorem 12.11 in which differentiability with respect to vector space 
topologies is treated. 

Theorem 7.2. Let (lUty and (Lj^) be satisfied and assume that g : R n — > R is bounded and 
differentiable. Moreover, suppose that j^(t,x,y,z) is Lipschitz continuous in x. Then for all 

x 1 — > (X x , Y x , Z x ), is differentiable, and 



p>2, the function R n -> 7^(R n ) x ^(R 1 ) x U>( 
the derivative is a solution of the BSDE LS2\) . 

Proof. By standard results, the mapping R n — > ^(R 1 ), x 1— > X x has a continuous version 
(which we assume being identical to the given one), and for all p > 1 there exists a constant 



C G R+ such that for x, x' G R n we have E(\X X -X?\ 2 p) < C 



x — x 



r\2p 



. See for example Lemma 



4.5.4 and Lemma 4.5.6 in [Kun90]. In order to be able to apply Theorem 12.11 we need to verify 
Condition (QT]). For this purpose, note that 

2 P -\ 



E 



\l(t,Xf,Y t x ,Z, 



Kt,X?,Y t x , 



Zf)\dt 



< E 



L\X X 



X? \dt 



2p 



< C(L,T, C)\x 



x 



l\2p 
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where L is a Lipschitz constant of I. Similarly, 

E[IJ \^-l(t,X x ,Y x ,Z x )-^-l(t,X x ',Y x ,Z x )dt\j ] < C'\x-x'\ 2p . 

This proves (CH]). Moreover, notice that = g{X?f) satisfies Condition (CG]). Thus the 

statement follows from Theorem 12.11 □ 

If in addition Condition (C(3]) is satisfied, we again obtain a sharper result stating pathwise 
continuous differentiability of an appropriate modification of the solution process. 

Theorem 7.3. Assume that (LU$, (T%M and (ZHJ), and suppose that the partial derivatives of I 
in the variables (x, y, z) are Lipschitz continuous. Then there exists a function £1 x [0, T] x R n — > 
(oj,t,x) i — > (X x ,Y x ,Z x )(u), such that for almost alluj, X x and Y x are continuous in 
t and continuously differentiate in x, and for all x, {X x ,Y X ,Z X ) is a solution of h31\). 



Let M > be a constant such that g, the derivatives of g, b and a, and the partial derivatives 
of I in (x,y,z) are all bounded by M. For all x G K let (X x ,Y X , Z x ) be the solution of the 
FBSDE (|31|) . To correspond formally to Theorem 12. 2\ in the setting of our FBSDE we have to 
work with 

l(u,t,x,y,z) = l(u,t,X?(u),y,z). 

But this functional fails to be globally Lipschitz in x. This is why we have to modify slightly 
the proof of Theorem 12.21 an d cannot just quote it. We start by showing that £(x) = g(X^) 
satisfies Condition (C3). 



Lemma 7.4. For all x G R n , h and i G {l,...,n}, letC Xi = \{g(X x+hei ) - g(X% )). Th 
for every p > 1 there exists a C > 0, dependent only on p and M , such that for all x, x' G 1 
and h, h! ^ ; 



E 

Moreover, for all t G [0, T], 

E 



|£*,A,i _ £x',A',i|2p < C Q X _ + \ h _ h 



'I 2 _L \h _ h'\ 2 W 



\C' h '' 1 - g\X?£ )^-X x T M < C{\x - x'\ 2 ) p . 



Proof Note that by Ito's formula g(Xf) = g(X§) + f* Vg(X x ) ■ a(s,X x )dW s + f* Cgds. Thus 
g(X x ) is a diffusion with coefficients a(s,x) = Vg(x) ■ a(s,x) and b(s,x) = Y2i=i ^i{ s i x ) %^ + 

\ X^j=i a ij( s i x ) dx 9< dx ■ ^ fiW) we have a G B lxd and b G B lxl . Therefore, by using standard 
results on stochastic flows (see Lemma 4.6.3 in |Kun90| ). we obtain the result. □ 



Proof of Theorem \7.3[ First note that it is well-known that X x may be chosen to be continuous 
in t and continuously differentiable in x (see for example Theorem 39, Ch. V, [Pro04j ). In order 
to prove that Y x has such a modification as well, note that Lemma f7.4l implies that £(x) = g{X?f) 
satisfies Condition (C3). Now let again U x,h = i(Y x+hei - Y x ), V x ' h = ^(Z x+he ^ - Z'f) for all 
x G M. n and h + 0. If h = 0, then define U x '° = £-Y x , V x '° = £-Z x and ( x '° = £-g(X£). It 
is enough to show that for all open bounded sets O C M. n+1 there exists a constant C such that 
for all (x, h) G O 



E 



sup \ u ^_ u x',h'^ P < C {\x-^\ 2 + \h-h'\ 2 y. 

te[o,T] J 
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As in Lemma 16.31 we will derive this estimate from Inequality (|30p . Notice that the assumptions 
of Theorem 17.31 guarantee that all the terms appearing in (|30f) . satisfy the same properties 
and thus provide the same estimates. There is one essential difference which is due to the 
appearance of X x instead of x in the first component of the line described by the integral 

lx,h(P) = {X x + 9(X x+he < - X x ),Y t x + 9(Y t x+hei - Y t x ), Zf + 6(Z x+hei - Z'f)). We therefore 
have to replace the prior A x t ' h by A x t > h = B^ h \{X x + he ^ - X x ) with B x t > h = ^<a x , h {0))d9 . 
Therefore we only need to show for (x, h) € O 

E[( [ T \Af h -A x '> h '\ds) 2pq2 ]^ <c(\x-x'\ 2 + \h-hf ) p . (33) 
Jo 

In fact, with A x > h = \(X x+he * -X x )\ih^Q and A x >° = VP, we have 

E( [ T \Af h -Af h '\ds) 2pg2 
Jo 

< / \A*> h \\B? h -B?> h '\ds+ / \Af h - A x '> h '\\B x '> h '\ds. 
Jo Jo 

The first summand satisfies 
rT 

.x,h\\r>x,h r>x',h' I i ^ 2p<? 2 



E( / \A x ' h \\B x ' h - B x '' h '\dsf 
Jo 

< E( jf \Af h \ J 1 |g( 7 ^W) - ^{ lx , >h ,{9))\d9dsf^ 

< (e( £ \Af h \ds) Apq2 y (e(£ J\i x m-^'^)ms)^ 

Lemma 4.6.3 in (Kun90j implies sup^ t h)eo ^( Jo" \As' h \ds) 4p9 < oo. Besides, 



T r l 



E( / / \l x , h (e))-j xl>hl {9)\d9ds) 



4pg 2 



JO 

<CE( sup \Xf' -X x \ + sup \xf +h ' ei - X x+he ^\ + sup \Y? - Y x \ 
v te[o,T] te[o,T] te[o,T] 

+ sup \Y t x ' +h ' ei -Y t x+hei \+ [ T (\Z X ' - Z x \ + \Z x ' +h ' e ^ - Z x+he '\)dtY Pq2 . 

t£[0,T] JO ' 

From this we can easily deduce E [(J^ \A X S > h \ \B X ' h - B x ' > h ' \ds) 2pq2 } < C{\x - x'\ 2 + \h - h'\ 2 )P. 

Similarly, E [( J Q T \A x ' h - A x ' > h ' \\B X ' ' h ' \ds) 2pq2 ] < C(\x - x'\ 2 + \h- h'\ 2 )P, hence © follows 
and the proof is finished. □ 



8 Malliavin differentiability of quadratic BSDEs 

In this section we shall ask for a different type of smoothness for solutions of quadratic BSDEs, 
namely differentiability in the variational sense or in the sense of Malliavin's calculus. Of course, 
this will imply smoothness of the terminal condition in the same sense. If the terminal condition 
is given by a smooth function of the terminal value of a forward equation, it will also involve 
variational smoothness of the forward equation. 
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Let us first review some basic facts about Malliavin calculus. We refer the reader to |Nua95j 
for a thorough treatment of the theory and to |KPQ97| for results related to BSDEs. To begin 
with, let C£ D (W nxd ) denote the set of functions with partial derivatives of all orders defined on 
M. nxd whose partial derivatives are bounded. 

Let S denote the space of random variables £ of the form 



Z = F((f Q h¥dWl) x <i< n , ■■■,(! hfdWt) X <i< n ) 



where F G C£°(R nxd ), h 1 , ■ ■ ■ ,h n G L 2 ([0,T];R d ). To simplify the notation assume that all h? 
are written as row vectors. 

If £ G S of the above form, we define the <i-dimensional operator D = (D 1 , ■ ■ ■ ,D d ) : S — » 
L 2 (Q x [0,T]) d by 



Dfc = J2d—{ hjdWt,..., KdW t )h^, 0<9<T, l<i<d. 
j—l ** JO 

For ( G 5 and p > 1, we define the norm 

neii 1;P = (E[i^ + (^ T |^ei 2 d0) i " 1 



It can be shown (see for example |Nua95] ) that the operator D has a closed extension to the 
space B 1,p , the closure of S with respect to the norm ||-|[i jP . Observe that if £ is ^—measurable 
then Dq£ = for 6 G (t, T] . 

We shall also consider n-dimensional processes depending on a time variable. We define the 
space hf p (R n ) to be the set of R n — valued progressively measurable processes u(t, u)te[o,T],wen 
such that 

i) For a.a. t € [0,T], u(t, •) G (B 1,p ) n ; 

ii) (t, u) — > DQu(t,uj) G (L 2 ([0,T])) dxn admits a progressively measurable version; 

\n\i, P =n( ; T Kt)i 2 dt) f + ( / T r T |z^iW) § ]i <^ 



Here, for y G W dxn we use the norm \y\ 2 = ^ jiUij) 2 - 
We also consider the space 

ro 1 ' 00 = n p> iD 1 ' p . 

We cite for completeness a result from [Nua95j that we will use in the next section. 

Lemma 8.1 (Lemma 1.2.3 in [Nua95j). Let {F n , n > 1} be a sequence of random variables in 
D 1,2 that converges to F in L 2 (Q) and such that 

sup E[\\DF n \\ L 2] < oo. 

nGN* 

Then F belongs to D 1 ' 2 , and the sequence of derivatives {DF n , n > 1} converges to DF in the 
weak topology of L 2 (£l x [0,T]). 

Let us now consider the BSDE 

Y t = Z- [ T Z s dW s + [ T f(s, Y s , Z s )ds. (34) 
Jt Jt 

Our assumptions on driver and terminal condition this time amount to 
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(El) / : 17 x [0, T] x K x M. d —* M is an adapted measurable function such that f(u,t,y,z) = 
l(u,t,y,z) +a\z\ 2 , where l(ui,t,y,z) is globally Lipschitz and continuously differentiable 
in (y,z); for all p > 1 we have E P [(/ Q T |Z(w, t, 0, 0)| 2 ds) 2 P] < oo; 

(E2) for all (t, y, z), the mapping $7 — > M, w i— > Z(u, i, y, z) is Malliavin differentiable and belongs 
to Lf p (R) for all p > 1. 

For any (u,t,y,z) and 6 € [OjT], the (a.e. valid) inequality holds true 

|D 9 Z(w,t,i/^))| < ^(w,t) + K e (u,t)(\y\ + |z|) 
where .Kg and Kg are positive adapted processes satisfying for all p > 1 



E[( sup / \K 6 (t,uj)\ 2 d0) P ) < oo and E[f / / |/? e (t, w)| 2 d0dt) P ] < oo 
v te[o,T]Jo 7 Wo 7o y 

Ivl.OO 



(E3) the random variable £ is bounded and belongs to 

We first consider the case where the terminal variable has no further structural properties, 
such as depending on the terminal value of a forward equation. For notational simplicity we 
shall treat the case of one dimensional z and Wiener process and so may omit the superscript i 
in D l etc. We will this time use the typical Sobolev space approach, hidden in Lemma 18.11 to 
describe Malliavin derivatives, which are in fact derivatives in the distributional sense. In this 
approach we shall employ an approximation of the driver of our BSDE by a sequence of globally 
Lipschitz continuous ones, for which the properties we want to derive are known. 

Let us therefore introduce a family of truncated functions starting with describing their 
derivatives by 

{— 2n , z < — n 
2z ,\z\<n 
2n , z > n. 

Then we have g n (z) = z 2 for \z\ < n, g n (z) = 2n\z\ — n 2 for \z\ > n, and thus |g r n ( 2 )| £ z 2 an d 
g n {z) — > z 2 locally uniformly on R for n — > oo. A similar statement holds for the derivative of 
g n (z): |<7^(.z)| — 2\ z \ an d g' n { z ) ~^ 2z locally uniformly on R for n — ► oo. 

With these truncation functions we obtain the following family of BSDEs: 

Y t n = e - [ T Zyw s + [ T [l(s, Y s n , Z n s ) + ag n (Z?)]da, n G N*. (35) 
Jt Jt 

From Proposition 2.4 of [KobOO] we obtain that there exists (Y S ,Z S ) G K°°(R) x L 2 (M) such 
that Y s n -» Y s uniformly in [0,T] and -» Z s in L 2 (M). 

Since the truncated equations have Lipschitz continuous drivers, Proposition 5.3 of |KPQ97| 
guarantees that (Y t n , Z™) G B 1,2 x D 1,2 with the following Malliavin derivative 

D e Y t n = and DgZ? = 0, if t G [0, 9), 

D e Y t n = D e i + ^ [d y l(Y?,ZZ)DgY? + d z l(Y s n ,ZZ)DgZ2 

+D e l(s, y s n , Z n s ) + ag' n (Z n s )DgZ n ^ ds - J D e Z?dW s , if t G [9, T\. (36) 

Now we aim at showing that the sequences DY n and DZ n are bounded in B 1 ' 2 , in order to 
use Lemma 18.11 This will be done by deriving a priori estimates in the style of the preceding 
sections, this time uniform in n. We therefore first show boundedness relative to the auxiliary 
measures Q n :=£(aj g' n {Z n )dW) ■ P, in the form of the following a priori inequality. 
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Lemma 8.2. Let p > 1. If the driver and terminal condition satisfy hypotheses o,nd 
(LWi), then the following inequality holds for the BSDE [36]) : 



E 



Qn 



sup / \D e Y t n \ 2 de) + E Q " 
te[o,T] Jo ' - 1 

■T 



T r-T 



D e Z r s l \ 2 d9ds 



< CE Qn 



o Jo 

T ,-T 



JO 



\D e l(s,Y s n ,Z2)\ 2 d9ds 



Proof. We will derive these estimates by proceeding in the same fashion as for Lemma 13.21 
Again, C%, C2, . . . are constants depending on the coefficients and p. 

Applying Ito's formula to e' 3 * \DoY t n \ 2 , using Equation (f36|) and simplifying as we did in the 
former sections we obtain (choosing (3 = M 2 + 2M) 



e 0t \D Y t n \ 2 + j\P s (A4\D e Y s n \ - \D e Z™\) 2 ds 



<e pT \Dotf -2 I e^D e Y s n D d Z?dW s + 2 I e ps \D e Y^D e l{s,Y^ Z?)\ds, (37) 



T 



where Wt = Wt — J * ag' n {Z™)ds, t £ [0, T], is a Q n — Brownian motion. 

We remark that since Esup te[0jT] \D e Y t n \ 2 < 00 we have E (f^ \D Y s n \ 2 \D e Z™\ 2 ds)^ 

and hence the process // eP 8 D g Y?D e Z?dW a is well defined. 
From (f3"T|) we obtain by taking conditional Q n -expectations 



< 00 



\D e Y t n \ 2 < CE^ \D e i\ 2 + \D B Y?\ \D e l(s, Y s n , Z^)\ds \T t 
L Jo 

Next, integrating in 9, using Fubini's Theorem and Doob's L p — inequality, we get 

r-T 



< CE^ n 



sup 

t<=[0,T] JO 



\D e Y t n \ 2 d9 



< CE Qn 



sup (E u 

t£[0,T] 
T 



T 



L J0 



T ,-T 



\Dotvm + 







\D S\ 2 dO + 

•T ,-T 
'0 JO 



JO 



\D e Y s n \\D e l(s,Y s n , Z^)\d9ds\T t 



\D e Y s n \\D e l( S ,Y s n ,Z^)\d9ds 



The last term on the right hand side of the preceding inequality can be simplified using Holder's 
and Young's inequalities with the result 



T r T 



\D e Y?\ \D 9 l(s, Y s n , Z^)\d9ds 



Jo 



< 



T _ , r T 



\D B Y Q n \ 2 d9 



I , fT 
2 



D e l(s,Y s n ,Z™)\ 2 d9 



ds 



< sup 
te[o,T] v Jo 



T 



\D e Y t n \ 2 d9 



I rT , rT 

2 



JO 



D e l(s,Y s n ,ZZ)\ 2 d9)'ds 

T ,-T 



<-^sup f \D e Y t n \ 2 dO + C 2 f f \Del(s,Y s n ,Z?)\ 2 d9d 
u i te[o,T] Jo Jo Jo 
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Since for a, b > we have (o + b) p < C^(a p + IP), by choosing C\ conveniently we obtain 



sup f \D g Y t n \ 2 d6 
te[o,T] Jo 

T 



< C&P" [( J' \D e t\ 2 &e) P + (f £ W(s, Y s n , Z^dsdeJ' 



(38) 



which provides the desired bound for the part of the Malliavin derivatives of Y n . Concerning 
the inequality for the Malliavin derivatives of the Z n part, we consult again Equation ([37]), from 
which we derive 



eP s \DgZ?\ 2 ds < e 



e^DeY^DgZ^dWs 



m\ 2 -2 [ 

Jo 

+2 / e l3s \D Y s n \\Dgl(s,Y s n ,Z^)\ds + 2 / Me Ps \D e Y™\\D e Z™\ds. 
Jo Jo 



Further estimate 



2 I Me^DgY^DgZ^ds < 4M 2 [ e? s \D 9 Y?\ 2 ds + J / e Ps \D e Z^\ 2 ds, 
Jo Jo 2 j 

f T f T f T 

2 / e Ps \D e Y s n \\D l(s,Y s n ,Z2)\ds < / e^ s \D e Y s n \ 2 ds + / e?'\D g l(t, Y s n , Z™)\ 2 ds. 
Jo Jo Jo 



Hence the initial estimate leads to 



\f e^DgZ^ds < e l3T \D e Z\ 2 -2 [ e^D Y^D e Z^dW s 
1 Jo Jo 

rT fT 

+(1 + 4M 2 ) / e f3s \D e Y s n \ 2 ds + / e 0s \Dgl(s, Y s n , Z™)\ 2 ds. 
Jo Jo 

Now for p > 1 integrate in 9, take Q n — expectations, using Fubini's Theorem as well as a 
stochastic version of it to estimate 



E y " 



T r-T 

Jo 



D e Z™\ z d9ds 

T 



< 



+ 



P 



+ ( sup 

-te[o,T] Jo 



\D e Y t n \ 2 d9 



T i-T 



JO 



D e l(s,Y s n ,Z2)\ 2 d9ds) + 



T i-T 



JO 



DgY s n DgZ n s d9dW t 



}• (39) 



We estimate the last term using Burkholder-Davis-Gundy's inequality, which results in 



E 



Qn 



T ,-T 



JO 



D e Y s n D e Z™d9dW i 



< C fi E 



Qn 



^ [j^ \D e Y?\\D e Z™\d9} 2 ds) 1 



Using Cauchy-Schwarz' inequality, we estimate further by 



T 



T 



\D e Y?\\DgZ™\de < ( / \D e Y?\ 2 d9 



\D B Z n s \ 2 d9 
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Then, with another application of Young's inequality, we obtain 

cT 

' s 1 1 o^s 







T ^ 2 

D Y?\\D Z?\dej ds 



< E c 



< E 



T 



T l 

\D e Y?\ 2 dey 



lA^lfdfl} 



ds 



{ sup f T \D e Y t n \ 2 de\ l { f T f T ' \D Z?\ 2 deds}* 

1 *e[o,T] j o J Wo jo J 



1 



c 7 



{ 



sup 

te[o,r] Jo 



T 



|D Y t n | 2 d0 



+ C 7 E Qn 



T 



JO 



iDsZ^dOds] 



Using this last estimate with C 7 chosen properly in conjunction with (|38[) in (|39p . we obtain 



r 

Jo 



|Dj^| 2 deds 



<C 8 E«"[(^ \D e Z\dey+(^ J \D e l(s,Y s n ,Z's)\ 2 ded s y]. (40) 

Combining inequalities (I38p and (I40p yields the desired estimate and proves the Lemma. □ 

In the same fashion as in Section [3j we can now combine the result of the a priori inequality 
under Q n ,n £ N*, with the inverse Holder inequality in disguise of Lemma 11.21 to upgrade the 
a priori estimates to the following one. In fact, we observe |<7^(z)| < 2\z\ for z £ R, n £ N*. 
Moreover, a careful analysis of the demonstration of Lemma 1 of [Mor07j shows that each 
J Z n dW is also BMO and there exists a constant K such that 



sup || I g' n (Z n )dW\\ B M0 2 < sup || f Z n dW\\ BM o 2 + \\ [ ZdW\ 

n£N* J n€N* J J 



BM0 2 ~- 



K < 00. 



(41) 



So by Lemma Fl.21 there exists a 1 < r such that £(J ag' n (Z n )dW)T and £(J aZdW)x are in 
L r (P) for all n € N* with r-norms bounded in n. So, again we may apply the argument based 
on the third statement of Lemma 11.21 uniformly in n. This, together with a similar argument 
applied to the sequence (Y n , Z n ) leads to the following a priori estimate. 

Lemma 8.3. Let p > 1 and r > 1 such that £(J ag' n {Z n )dW)T £ L r (P) for all n £ N* with a 
uniform bound. Then there exists a constant C > 0, depending only on p, T and K (from \41\ ), 
such that with the conjugate exponent q of r we have 



E" 



■T 

JO 

< c{e p 



\D e Y t n \ 2 d6dt 

T 



+ E i 



T ,-T 



\D e Z™\ 2 d6ds 







Jo 

T r T 



JO 



K e (uj,t)\ 2 d9dt) 



+E } 



sup 

te[o,T] Jo 



\K ft (ojA)\ 2 d9 



|e| 2 + / 



Ms 



\ 4p<j 



4, 1 
2? 



< OO. 



Proof. The proof of the lemma is achieved in three steps. 

We start by applying the third statement of Lemma 11.21 as in the proof of Theorem 13.11 to 
the result of Lemma l8.2| from which we obtain with a constant C\ not depending on n 



E 1 



T r-T 
JO 

< dE p 



D e Y t n \ 2 d9dt 

T 



+ E i 



T rT 







T 

Jo 



Jo 

T 



\D e Z™\ 2 d6ds 



\D e l(s,Y s n ,Z™)\ 2 d6ds 



i»r 
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In a second step, we have to estimate the last term of the preceding equation. From Condition 
(E2) we obtain with another universal constant 



T i-T 

\D e l(s,Y s n ,Z^)\ 2 d0dt 



JO 



<C 2 [ T f T \[K e {u;,t)\ 2 + \K e {u:,t)\\\Y t n \ 2 + \Z^\ 2 )]d6dt 
Jo Jo 

<C 2 { / ( \K e (uj,t)\ 2 d8dt+ sup f \K e (u,t)\ 2 d0 [ \Y s n \ 2 + \Z^\ 2 ds\. 
L Jo Jo te\o,T] Jo Jo J 



te[o,T] Jo 

Hence an application of Holder's inequality results in 

cT r-T . „„2 , r-T r-T . „„2 



E[(jT jf |D e /( s ,y™,Z s m )| 2 d0d S ) P9 ] <C 3 E[(^ jT |K (o;,i)| 2 d0d S ) M 

+C 4 E[( sup / T |^(a;,t)| 2 d^ 2M2 ]5E[( /^/f + |Z^W) W ]i (42) 
v ie[o,TlJo y Wo y 



(To \ 
Jo l-^TI + ds] ] term. For 

this purpose, we shall use another application of Theorem 13.11 uniformly in n. It requires 

the intervention of a different family of measure changes depending on n, which can again be 

controlled by the BMO property of the intervening martingales and the third statement of 

Lemma 11.21 In fact, comparing (j35[) with ([5]), we see that the analogue of H has to be given by 

h n (Z n ), where 

g n (z) \ z ,if \z\ < n 



h n (z) :-- 



2n\z\—n 2 ■ r i i _ 
— — , if \z\ > n, 



which is obviously well defined for all z. In this situation, the stochastic integrals of h n {Z n ) 
generate BMO martingales with uniformly bounded norms. More precisely, since g n (z) < z 2 , 
we have sup ngN * |/i n (^)| < \z\, z £ R. A careful analysis of the demonstration of Lemma 1 of 
[Mor07] shows that each Z n is also BMO and there exists a constant K such that 

sup || / Z n dW\\ B Mo 2 + \\ / ZdW\\ B Mo 2 = K <oo. 

nGN* J J 

Due to the definition of h n , we may extend (I41D to 

sup || / h n (Z n )dW\\ BM o 2 < sup || / Z n dW\\ B Mo 2 < K. 

nSN* J nGN* J 

By Lemma 11.21 we may assume that £ (J ah n (Z n )dW)T and £(f aZdW)x are in L r (P) for 
all n £ N* with r-norms bounded in n. So, again we may apply the argument based on the 
third statement of Lemma 11.21 uniformly in n. We obtain from Theorem 13.11 with the settings 
A t = J(w,t,0,0), l(u,t,Y t n ,Z?) := 1(oj, t, Y t n , Z™) - l(u, t, 0, 0), ( = £ and H t = h n (Z?), for all 
7 > 1 the inequality 

su P E p [( sup |y t n | 2 + [ \Z?\ 2 dsY] < C 4 E P [(|£| 2 + [ T \l{oj,s,0,0)\ 2 dsY q ]** . (43) 

nGN* L M6[0,T] JO ' - 1 L\ Jo ' - 1 

Plugging this inequality into (|42l) terminates the proof of the Lemma. □ 
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Our main result can now be proved. 



Theorem 8.4. Assume that driver and terminal condition satisfy hypotheses (I^jty, (7Hj) and 
Then the solution processes (Y t ,Z t ), of (3$ belongs to B 1 ' 2 x (B 1>2 ) d and a version of 
(DgYf, DgZt) satisfies for 1 < i < d 



D l e Y t 
DlYt 



0. 



D\Z t = 0, 

T 



m+ 

i-T 



*e [0,0), 

d y l(s, Y s , Z S )D\Y S + 8 z l(s, Y s ,Z s )D g Z s + D^s, Y s , Z s ) + 2aZ s D e Z & 

D l e Z 8 dW 8 , te[6,T}. 
Moreover, {DtYt : < t < T} is a version of {Zt : < t < T}. 



ds 
(44) 



Proof. Again, we simplify notation by just considering the case of one dimensional Z and 
Wiener process. We first apply Lemma 18.31 to obtain a bound of the L 2 -norms of the pro- 
cesses (DgY™ , DgZ™), which is uniform in n. Lemma 18.31 furthermore allows an appeal to a 
weak compactness result to deduce the existence of a pair of processes (Uo,t, Vo,t), < 9,t < T, 
and a subsequence (n«) such that (DgY^ li (u), DgZ^ i (u)) converges to (Ug tt (uj), Vgj{w)) m the 
weak topology of the L 2 space of random variables with values in L 2 ([0, T] x [0, T]). 

For almost all t, Lemma 18 . 1 1 implies that (Yj, Zt) is Malliavin differentiable and the equality 
(DgY t , DgZ t ) = (Ugj, Vgj) holds almost everywhere in Q x [0, T]. 

It remains to use these convergence properties to deduce convergence term by term in ([36]) 
to (PI). 

We first show that the stochastic integral terms converge weakly in L 2 (Q). To this end 
let ^ £ L 2 (f2) be .^-measurable. Then there exists a predictable ip G L 2 (!J x [0, T]) with 
* = E(^) + | T ip s dW s , and hence 















lim E 


^ f T D g Z^dW s 


= lim E 


[ 7 \/; s DgZ^ds 


= E 


/ i) s D e Z s ds 




. Jo 




Jo 




Jo 



E 



/ D e Z s dW s 



which shows that J DgZ^dW s converges weakly to J D e Z s dW s in L 2 (ft x [0,T]). 
Next observe that for any bounded ^r- m easurable random variable B we have 

E[B [ T (g'(Z^)D e Z^ - g>(Z t )D Z t )dt] 



B f \g' ni (Zr) ~ g'(Z t ))D d Z^dt]+E[B C g'(Z t )(D e Z^ - DgZ t )dt] (45) 
Jo Jo 



E 



The first summand on the RHS of Equation (|45p is bounded by 



(esssup \B\) sup ( E / D e (Z^) 2 dt 



E 



(g'(Zn-g'(Z t )) 2 dt] 



which converges to as i — > oo. The second summand on the RHS of Equation (|45f) converges 
also to since Bg'(Zt) € L 2 (0 x [0, T]) and DgZ™ 1 converges weakly to DgZf. Since B was 
arbitrary we have shown that Jq ^.{Z^^DgZ^dt converges to g'(Zt)DgZtdt in the weak 
topology of L 1 ^ x [0,7]). 
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Finally we come to the various derivative terms of I. The boundedness of the partial deriva- 
tives of I as well as Condition (Ej2|) and Inequality (|43f) imply 
cT f-T 



supE / / [D e (l(s,Y s n ,Z^)] 2 dsdd 

n JO JO 

< dsupE [ T C ' \d y l(Y^,Z2)D e Y^\ 2 + \d z l(Y s n ,Z^DgZ n \ 2 + [(Del)(s,Y s n ,Z^} 2 dsd9 

n Jo Jo 

< C 2 {supE [ [ \D e Y s n \ 2 + \D e Z^\ 2 + \K e (s)\ 2 dsd0 

•> n Jo JO 

sup \K e (s)\ 2 d6 [ T (\Y s n \ + \Z?\) 2 ds } 
efo,Tl Jo J 



+ supE 



< oo. 



se[Q,T] 



Thus, by LemmaEEl for almost all s G [0,T], l(s,Y s ,Z s ) belongs to B 1,2 , and D e (l(s, YJ 1 , Z?)) 
converges to D e (l(s,Y s ,Z s )) weakly in L 2 (Fl x [0,T]). Since the partial derivatives of I are 
continuous we have D 8 {l{s, Y s , Z s )) = d y l(Y s , Z s )D e Y s + d z l{Y s , Z s )D e Z s + (D l)(s, Y s , Z s ). □ 

We next assume more structural properties for the terminal variable. More precisely, we will 
turn to the framework of forward- backward systems. Given a d— dimensional Brownian motion 
W and an x £ R ra , we denote by Xt = {X} , • • • , X™) the forward part given by 

X t = x+ [ b(s,X s )ds+ f [a(s,X s )]*dW s . (46) 
jo Jo 

The coefficients are supposed to satisfy 

(PI) b, Gi : [0,T] x W 1 — > W 1 , are uniformly Lipschitz; b(-,0) and o"j(-,0) are bounded for 
1 < i < d; a(t,Xt) is a d x n matrix and [cr(t, Xt)]* represents its transpose. 

For the backward part we consider 

Y t = g(X T )- [ Z s dW s + [ [l(s,X s ,Y s ,Z s ) + a\Z s \ 2 ]ds, t € [0, T] (47) 
Jt Jt 

where the driver and the terminal conditions are supposed to satisfy the following assumptions 

(P2) I : [0, T] x M. n x R x R d — > R is globally Lipschitz and continuously differentiable in (x, y, z); 
/(•, 0, 0, 0) and d x l(t, x, y, z) are bounded by a constant M. 

(P3) g : M. n — > M. is a bounded differentiable function with bounded first partial derivatives. 

Before stating the main theorem we recall that the processes X, Y and Z all depend on the 
variable x. For ease of notation we omit the corresponding superscripts. In this setting our main 
result is the following. 

Theorem 8.5. Suppose that the coefficients of the SDE J^ffi ) and the driver and terminal con- 
dition of the BSDE fl^7| ) satisfy conditions (PI), (P2) and (P3). Then the solution processes 
(X, Y, Z) possess the following properties. 

• For any < t <T, x G R, (Y t , Z t ) £ B 1 - 2 x (B 1 * 2 )^ and a version of {(D^Yt, B\Z t )\ < 
G,t < T} satisfies for 1 < i < d 

D\Y t = 0, DlZ t = 0, t€[0,9), 

D l e Y t = d x g{X T )DiX T + J [d x l(s,Y s ,Z s )D l e X s + d y l(Y s ,Z s )DiY s 

+d z l(Y s ,Z s )D l e Z s + 2aZ s D i e Z s jds- j\D l e Z s ,dW s ), te[9,T\. (48) 
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Moreover, {DfYt; < t < T} defined by the above equation is a version of {Zt; < t < T}. 

• The following set of equations holds for any < 9 < t < T and x G W 1 , P-almost surely, 

D e X t = d x X t (d x X e )- 1 a(9,X e ) 
D e Y t = d x Y t (d x X e )- l a{e,X e ) 
Z t = d x Y t (d x X t )- l a(s,X t ); 

and DgZt = d x Zt(d x Xg)~ 1 a(9, Xg) for almost all (u,t). 

Proof. Theorem 2.2.1 of [Nua95j assures existence, uniqueness and Malliavin differentiability 
of solutions of SDE (|46|) under Hypothesis (PI). Moreover the solution processes satisfy Xt G 
(D 1,00 ) n for any t G [0, T] and 1 < i < d and the following equation holds: 

D\x t = o, t€[o,e), 

D'gXt = <r(9,X g )+ I d x b(s,X s )D i g X s ds+ [ d x a{s,X s )DiX s dW s , te[e,T], 
Je Je 

(see f.ex. Theorem 2.2.1 of |Nua95j ). 

Let us next check the validity of hypotheses (E1)-(E3) for the driver of our BSDE, for 
simplicity in the one dimensional case. From Condition (P3) it follows that the function g and 
its derivative are bounded. In combination with the fact that Xt G D 1 ' 00 this implies that 
g(Xx) G B 1 ' 00 , i.e. in the setting of Theorem 18.41 Condition (E3) is verified. From Condition 
(P2), we have \l(t,X t ,0,0)\ < M(l + \X t \) for t G [0,T]. The fact that X t G B 1 ' 00 then entails 



E[( J \l(t,X t ,0,0)\ 2 ds) ] < oo for all p > 1. Hence (El) is satisfied. Condition (P2) includes 



the statement that d x l(-,x, -, •) is bounded. Therefore we have 

\D e l(t,X t ,y,z)\ = \d x l(t,X t ,y,z)D e X t \ < M\D 6 X t \ 

with some constant M. Using the fact that X t G B 1 ' 00 we obtain J Q T \M D e X t | 2 d6»dt^] < 

oo for p > 1, which means Condition (E2) also holds. 

With conditions (El), (E2) and (E3) verified we can apply Theorem 18. 41 which implies the 
Malliavin differentiability of (Y S ,Z S ) and proves the first block of results. 

For the second part of the theorem, the representation formula of DX is standard (see f.ex. 
Chapter 2.3 in [Nua95] ) . The representation of Z by the trace of DY being granted, we only 
have to prove the representation formulas for DY and DZ. For this purpose, we apply Ito's 
formula to d x Y t (d x Xg)~ 1 a(6, Xg), then use (J1J) to represent the d x Y t term. We further use the 
representation of DX to account for the terminal condition. This way we obtain (|48p with 
DgY t = d x Yt(d x Xo)~ l o-(9, Xq) and DgZ t = d x Zt{d x Xg)~ 1 a{9,Xg). The representation follows 
from uniqueness of solutions for the BSDE. □ 

Example: 

We finally study a specific setting of Theorem 18.51 We assume that £ and D^ are bounded by 
M. Assume further that the driver does not depend on y. Then, choosing 9 = t, representation 
IBj) can be simplified to 



Z t = £>*£+ / \d z l{Z s )D t Z s + 2aZ s D t Z s ]ds- [ D t Z s dW s . 

Since j ZdW G BMO and if we further assume J d z l(Z)dW G BMO we may change the mea- 
sure to Q = Si f d z l(Z) + 2aZdWT) P- Hence we obtain, by applying conditional expectations 



Z t = E Q [D t $\T t ] < M¥P[l\F t } < M. 
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This means Z £ L°°. This way we recover the Malliavin differentiability results of [HM06J from 

our main result. 

Remarks: 

1. The methods of proof of this Section, building upon a truncated sequence of Lipschitz 
BSDEs, could also be used in the treatment of the differentiability problem in Section [6l This 
sequence would allow the use of the results in |KPQ97|, which, combined with the a priori 
estimates of sections [3] and 0] would imply differentiability. 

2. Our main results allow less restrictive hypotheses. For example in Section [3l we assume 
for our a priori estimations that ( € L p for all p > 1. An analysis of the proof clearly reveals 
that to obtain estimates in TV or L p we only need that for all p € (2, 2pq 2 ]. We chose to 
write ( € L p for all p > 1 not to produce an overload of technicalities in a technically already 
rather complex text. 
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